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Abstract 

We examine some excited state energies in the non-unitary integrable quantum field 
theory (IQFT) obtained from the perturbation of the minimal conformal field theory (CFT) 
model A^3,5 by its operator 02, i- Using the correspondence of this IQFT to the scaling limit 
of the dilute A2 lattice model (in a particular regime) we derive the functional equations 
for the QFT commuting transfer matrices. These functional equations can be transformed 
to a closed set of TBA-like integral equations which determine the excited state energies in 
the finite-size system. In particular, we explicitly construct these equations for the ground 
state and two lowest excited states. Numerical results for the associated energy gaps are 
compared with those obtained by the truncated conformal space approach (TCSA). 
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1 Introduction 



The problem of calculating the excited state energies in finite volume integrable quantum field 
theories (IQFT's) has recently received much attention In particular, in reference [|I| it was 

shown that the functional equations for the IQFT "commuting transfer-matrices" introduced 
in can be transformed into integral equations which generalise the standard ground state 
thermodynamic Bethe Ansatz (TBA) equations to excited states. The results of apply to 
massive IQFT's obtained from 0i 3 perturbations of the minimal conformal field theory (CFT) 
models, which are related to the Uq{sl{2)) quantum algebra. To include IQFT's related to 0i,2 

or 02,1 perturbations of these minimal CFT models, one has to generalise this approach to the 

(2) 

case of the g-deformed twisted Kac- Moody algebra A2 . Some results in this direction were 
obtained in |lTO|,p]T[]. 

In this paper we consider the IQFT obtained from the perturbation of the non-unitary 
(c = —3/5) minimal CFT A^3,5 by its operator 02, 1 with conformal dimensions A = A = 3/4. 
It is known that the particle spectrum of this model consists of three kinks of the same mass 
which interpolate between two degenerate vacuum states. The exact S-matrix of these kinks 



was conjectured by Smirnov [|I2[. A numerical study of the energy spectrum of this model 
was performed in |jl3| using the Truncated Conformal Space Approach (TCS A) fl^ . The TBA 
equations for the ground state energy were originally conjectured in |T^, and more recently, were 
derived in [16] from the exact S-matrix of the model. In this paper we derive the functional 
relation for the eigenvalues of the "QFT transfer matrices" which is, in principle, sufficient to 
determine the whole energy spectrum of the model. In doing this we do not directly follow the 
route of works but instead obtain the required functional relation as the scaling limit of 

the associated lattice model results. Fortunately, a suitable solvable lattice model whose scaling 
limit describes the above IQFT is known. It is a particular case of the off-critical dilute-y42 



lattice model [|18| in a certain regime. 

The organisation of the paper is as follows. In Section 2, we briefiy review the related results 
of where the TBA integral equations for the ground state energy of the model were derived 



using the conventional TBA approach [T^, The massive field theory limit of the related 
dilute-A2 lattice model is considered in Section 3. In Section 4, we present the functional 
relation for the QFT transfer matrices along with some numerical results on the positions of 
zeroes for the largest three eigenvalues of these transfer matrices. Using these results we then 
derive the TBA-like integral equations which correspond to the ground state and two lowest 
excited states energies of the finite-size QFT described above. Comparisons are made between 
the two lowest excited state energy gaps computed from the above TBA-like integral equations 
and corresponding results of the truncated conformal space approach (TCSA). 
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2 The /S-Matrix and Ground state TBA Equation 



Consider the massive QFT described by the action 



[2.1] 



where Amz 5 represents the action of the c = —3/5 minimal non-unitary CFT Al3,5, the operator 
02,1 has conformal dimensions (3/4,3/4), and (7 is a couphng constant of dimension (mass)^/^. 



It is known that this action defines an integrable QFT in the sense that it possesses an 
infinite number of non-trivial integrals of motion. The particle spectrum of the model consists of 
three kinks (of the same mass m) interpolating between two degenerate vacuum states (labelled 
and 1). The precise relationship between the perturbation parameter g and the kink mass m 
is O 

^r(^)r(|) 
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2V3nT (X) r (I) 



m2 



0.25300091957. 



.1171- 



The exact kink-kink ^-matrix (as first conjectured in |12]) can be written in the form [13 

R{9) W 



S'JM 




a,f3,-f,5e {0,1} 



[2.2) 



(2.3) 



where 9 is the physical rapidity variable, the parameters x = — ^, po = 2cos/x and pi = 1, 
and the normalisation factor is 



R{9) 



sin (f - f ) sin (f - f ) 



(2.4) 



sin + ^) sin (f + ^) ' 

Here, the functions W{...) denote the Boltzmann weights of the critical hard hexagon lattice 
model |2^ (in an unphysical regime). With a suitable normalisation these Boltzmann weights 
can be expressed as 

sin /i sin(2yU-|-'u) 



W\ 



w\ 



w\ 



u 



u 



u 



u 



u 



sin2/i sin(/i— -u) 




sm p 
sin 2p 



smw 



sin(/i— m) 



(2.5) 



sin/i sin(2/i— m) 
sin2/i sin(/i— -u) 

sin(yU-l-'u) 
sin(yU— -u) 
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where u = x9 is the lattice spectral parameter and /i = 37r/5 is the crossing parameter. 

The S'-matrix ( p.3|) defines a factorized scattering theory. This implies that a system of 
N kinks distributed along a large spatial circle of length L can be analysed via the usual 
Bethe Ansatz (BA) approach. The rapidities ^i, . . . , of the respective kinks are consequently 
constrained by the Bethe- Yang equations pO| 



^imL sinh 6k 



where A{9] 9i 



'Ai9k;9i,...,eN) = -l, k = l,...,N 
,9n) are the eigenvalues of the "scattering transfer matrix" 



N 



9 



N, 



{a} 



2 = 1 



CtiPi+1 



{9-9i). 



[2.6) 



[2.7) 



From equation ( p.3[ ) it is obvious that this is just the transfer matrix of the inhomogeneous 
lattice model with the Boltzmann weights given by ( |2.5| ). The eigenvalues of this lattice model 
were found in [^] using the analytic Bethe Ansatz approach. This result (together with some 
assumptions on the string structure of the associated Bethe- Ansatz equations) enables one to 
follow the standard TBA procedure |]T6| . Below we briefiy review these calculations. 

The states of the system in the thermodynamic limit N oo, L ^ oo, are specified 
by the densities of the rapidity distributions for the kinks and quasi-particles arising in the 
diagonalization of the transfer matrix 



As shown in |16], only one type of quasi-particle 
is required in our case. The densities of the rapidity distributions for the kink and the quasi- 
particle (and the associated densities of "holes" in these distributions) are then determined by 



the following integral equations |]T6 



m 
2^ 



where 



1 /"OO 

Sjfl coshe = aj{9) + aj{9) + ^ ^j,ki 

k=o 



o-k[ 



d9\ 



J = 0, 1. 



v^sinh(2^) 



(2. 



(2.9) 



vr sinh(36') 

Here cro{9) denotes the density of the kink rapidity distribution which is normalised as 

ao{9)d9=N, (2.10) 

1 

and ao{9) denotes the corresponding density of "holes" in this distribution. Similarly, <7i{9) and 
^i{9) denote the rapidity and hole densities for the quasi-particles. The density of states in the 
system is determined by the usual combinatorial entropy p5[ 



S[(ro{9),(ri 



fc=0 



[a,{9)+a,{9)] \n[a,{9)+a,{9)]-a,{9) \na,{9)-a,{9) \na,{9) \ d9. 

2.ir 
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The free energy of the system is then defined as the functional 

J^[aoie),a,ie)]=S[aoie)]-^S[aoie),aM (2-12) 
where 1/R acts as the temperature parameter, and the energy is 

/oo 

cosheao{e)de . (2.13) 
-oo 

Here, the parameter 

e = log(m/?o)/(87r) (2.14) 

is the vacuum energy (per unit length) for the QFT defined by the action ( p.l| ) |T5|,^, and Rq 
is a non-universal ultraviolet cutoff parameter. Using the standard calculations, one obtains the 
equilibrium free energy in the form 

/o(mi?) = - = e / cosh^ log 1 + e-'"^^^ dO (2.15) 

L 2-nR J_^ V y 

where the pseudo-energies 



8,(6) =\og [a,ie)/a,ie) 
are determined by the TEA integral equations [p!5| ,p!6[| 



ej{e) = SjflmR cosh{e) + J2 [ ^j,k{0 - O')^og (l + e'^'^^'A dO' , j = 0, 1. (2.16) 
As is well known, the free energy ( p.l5| ) can be re-interpreted as the ground state energy 

E,iR) = RfoimR) = "^""^ ^"^^"^"^"^ -^T coshO log f 1 + e'-^^) ) dO (2.17) 

of a finite-size system defined on a circle of circumference R. The leading asymptotics of Eq{R) 
in the ultraviolet limit i? — 0, can be calculated using the standard "dilogarithm trick" |p6|| . 
This results in 

E,(R) ~ (2.18) 

which is in agreement with the expected form Eq{R) ~ — (7r/6-R)(c— 24Ao), where Ao = — 1/20 is 
the lowest conformal dimension of the operator algebra associated with the minimal CFT model 

-^3,5- 



5 



3 The dilute- A2 lattice model 



The IQFT model considered above is associated with a certain regime of the solvable dilute-742 
lattice model [|1^,^. The general dilute-AA/ model is an interaction-round-a-face (IRF) model 



defined on the square lattice which has discrete spin variables (or "heights") located at the 
lattice sites. These heights can take one of M (with M > 2) possible integer values 1 < ctj < M, 
subject to the restriction that neighbouring sites either have the same height or differ by ±1. 
The explicit expressions for the Boltzmann weights of an elementary face of the lattice are given 
in Appendix A. They are parameterised through the elliptic theta functions ^j{u\ p) of (spectral) 
variable u and nome p = e~'^ which are also given in Appendix A. Moreover, these Boltzmann 
weights depend on an additional parameter A which should be chosen from the discrete set of 
values 

^ = MmWY ke{l,...,M, M+2, 2M+1}. (3.1) 

It is only when k = M ot k = M+2 that the model is physical in the sense that all of the 
Boltzmann weights are real and positive (for some range of the variable u on the real axis). 
Here we consider a specific unphysical case (corresponding to = 1 in (pTTP) with 



M = 2, p>0, 0<u<3\, A = 7r/12. (3.2) 

We show below that the scaling limit of this model is described by the IQFT considered in the 
previous sections {i.e. the IQFT obtained as a perturbation of the minimal non-unitary CFT 
TVIs^s by its operator 02, i)- 

The row-to-row transfer matrix (with periodic boundary conditions in the horizontal direc- 
tion) is defined by 

where {a} = {ai, 02, . . . , otv} and {b} = {61, 62, • • • , ^tv} denote heights on two consecutive 
rows of the lattice, periodic boundary conditions imply that Qn+i = ai and &Ar+i = bi, and 
is the number of sites per row. Since the Boltzmann weights of the model satisfy the Yang- 



u (3.3) 



Baxter equation |T8|, the transfer matrices T{u) with different values of the spectral parameter 



u commute among themselves (for fixed values of p and A) 

[T{u), Tiu')] = (3.4) 

and hence can be simultaneously diagonalized. The corresponding eigenvalues of the transfer 
matrix ( ^.31 ) are given by 



A{u) = UJ (^h{2X-u) h{3X-u)y 



Q{u-X) 
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+ {h{u)h{3X-u] 



N 



Q{u) Q{u-3X) 



+ iu~^ - h{u) h{X-u) 



Q{u-X) Q{u-2X)' 
N Q{u-AX) 



(3.5) 



Q{u-2X) 



where 



Q{u) 



N 



u 



(3.6) 



and where the function h{u) coincides (up to a simple factor) with the standard theta function 
defined in ( |A.1| ) 

h{u) = "p^^l^ i}\ {u]p) , p = e"'^ . 



The numbers {vj} are solutions of the following set of Bethe Ansatz (BA) equations 



UJ 



h{X+Wj] 
h{X-iVj] 



N 



N 

n 



h{ivj—ivk+2X) h{wj—Wk—X) 
^ h{ivj—Wk—2X) h{ivj—ivk+X) 



J 



N 



(3.7) 



(3.8) 



with = exp(z7rm/(M+l)), m = l, ...,M. 

Now consider the Hamiltonian of the associated one-dimensional quantum spin chain 



u=0 



+ const 



(3.9) 



where 5 is a parameter with the dimension of length. The corresponding energy eigenvalues 
(with a suitable choice of the constant term in (|3.9|)) can be expressed as 



E 



1 d 



N 



N 



LU- 



h{u + A 



IV j 



4:5 du \ h{u — A — ivj) 



(3.10) 



M=0 



The scaling limit of this lattice model can also be analysed within the TBA approach. When 
N —>■ oo the solutions {vj} of (|3.^ ) converge to certain asymptotic patterns in the complex v- 
plane, which can be viewed as collections of "strings" pS. An ^-string is a set of £ > 1 complex 
roots with the same real part, symmetric with respect to the real axis of v and equally spaced 
(with the spacing 2A) along the imaginary axis. The total number of roots £ in a string is called 
the length of the string. Numerical calculations suggest that for — oo the string spectrum 
of (|3.8| ) is saturated by the strings of lengths i = 1,2,3 only (in the sense that the number of 
strings of other types in a generic solution of (p.8|) grows slower than A^). Assuming this picture 
is correct, it is possible to write the roots vj which solve ( |3.8| ) in the form 



+ 2Af + 0(e 



1,2,3 



(3.11) 
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where 



A 



(3.12) 



and the (real) numbers o'^^ determine the centres of the strings. It is to be noted that for 
non-vanishing values of the elliptic nome p > the convergence parameter a in ( p.ll| ) is positive 
and non-zero. The integers N^^"* denote the total number of strings of length £ in a particular 
solution of (|3.8| ). As follows from (|3.6|) these numbers are restricted by 



^£ArW = N + o{N). 



(3.13) 



^=1 



In the thermodynamic limit N oo, the centres of the strings form continuous distributions 
and the BA equations ( p.8|) lead to the following integral equations^ for the densities of strings 
pe{e) and "holes" pe{e) [0 



ie) + psie) + J2 r B,,t{e-e')pt{9')d9\ s,t = i,...,3, (3.i4) 

t=l -^-2r 



where 



s 



fc=i 



s t 



Bs,m = E E h + ^^^^^ + '^-) - + + 

k=l m=l 

and where the function aj{9) is defined as 

1 d 



(3.15) 
(3.16) 



2m dO 



log 



h{X]/2 + iO/A) 

yh{\j/2-i6/A)\ • 



(3.17) 



Note that the functions aj{0) (as well as hs{0), Bs^t{d) and all of the densities p(,{9) and pe{0)) 
are periodic with period 4r. Also, the densities pi are normalised by the condition 



2t 



Pi{9)d9 = N^^^/N 



2t 



which from ( p. 13 ) implies that 



1. 



(3.18) 



(3.19) 



This approximation is valid provided Np^ ^ 1 in addition to ^ 1. In this section we will assume that 
this extra condition is also satisfied. 
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This relation together with equation (|3.14|) (with s = 3) imphes that P3{6) = 0, and hence it is 
possible to eliminate the density P3{6) from (|3.14|) . Re-denoting the remaining densities as 



where the parameter 5 is the same as in (|3.9|) , one obtains the integral equations 

1 /"OO 

5-H,flMo+'^r,T) = a,{e) + a,{e) + Y. / ^,^k{e-e\T) au{e') de\ j = o,i. (3.21) 

fc=o -^-"^ 

Here, the functions 

r) = (5,, - 5j,k-i - 5i,fc+i)0o(^, r) (3.22) 



5o 

n=— 00 



1 jg^.^ COsh(2AXn) -KTL 



4r cosh(6Ax„) 2r 

n=— 00 ^ ' 

and the parameter A = 7r/12. The largest eigenvalue Ao(m) of the transfer matrix (|3.3| ) (and 
hence the ground state energy of the Hamiltonian ( p.9| )) is determined by the solution of ( |3.8| ) 
that consists only of 1-strings. The corresponding density 

p'^r\d) = Md,r) (3.24) 

can be easily determined from ( ^3.211 ). Then from ( ^.51 ) one obtains 



, /^ r A / ^^ , /^^i (?/+A+z^/4)^9i (4A-n+ze/4) \ 

'^^") = J- [n J = ^^^^ + 1 ( V(A+.^/4)^,(4A+.^/4) ) 

(3.25) 

where 

p = h{2X)h{3X). (3.26) 
Finally, the spectrum of the Hamiltonian ( p.lO| ) for large N reads 



d 



46 du 



M = 



2r 



En = - — —\ogK{u) +2tiN 0o(^+2r, r) (Tq (0)^0 (3.27) 



-2t 



where the first term represents the vacuum energy and the second term represents the excitation 
energy. The gap in the excitation spectrum is then determined by the value of 0o(6'+2r, r) at 
^ = where the one-particle excitation energy has a minimum. 

Now consider the scaling limit. The leading singularity of (|3.25|) for p ^ is given by 



X2 sill 4?x 

log K(M)si„g = p'^ logp (3.28) 

vr 
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which determines the thermal exponent ot and correlation length exponent v to be 



2-a = 4, v = 1. (3.29) 

Let us now identify the parameter 8 in ( |3.9| ) with the lattice spacing constant, and take the 
limits N ^ oo and 6-^0, while keeping both the (dimensional) length of the chain L = N6 
and the correlation length Rc ~ p~'^6 finite. This means that 

p"^ ^ N~\ 6r^N-\ N^oo. (3.30) 

We shall also assume that L ^ (which is equivalent to Np'^ ^ 1, see footnote 1 on page 8). 
In the above limit, the gap m in the one-particle spectrum of (|3.27| ) {i.e. the mass gap) tends 
to a finite value. Indeed, one has 

27rrVo(^+2r,r) = mcosh^ + 0(/), P ^ 0, (3.31) 

where 



m 



aV3p^6-\ (3.32) 



and hence expression (|3.27| ) becomes 



^^^^^ m^Llog(mLo) ^^^ cosh^aoc^^. (3.33) 



61T 



It is to be noted that the bulk vacuum energy term here is determined entirely by the singular 
part (p.28|) of the largest eigenvalue of the transfer matrix. The parameter Lq is non- universal, 
and in particular, it depends on an overall shift of the energy spectrum in (|3.9|) . It is obvious 
that the energy expression (|3.33|) is identical to that given by ( p.l3|) and (|2.14| ) of the previous 
section. Similarly, using ( p.31| ) and 

M0,r)=m + Oip'), p-^0 (3.34) 

with 0(6*) given by (|2.9| ), it is easy to see that the integral equations (|3.21| ) in the limit ( p.30|) 
become identical to ( |2.(j| ). And finally, the density of states in this lattice model is, of course, 
determined by exactly the same combinatorial entropy functional as in (2.11|). 



4 Excited states energies for the finite-size system 

In the previous section we established an equivalence between the QFT (|2.1|) and the scaling 
limit of the dilute A2 lattice model. We have shown that the density of states in these two 
systems are described by the same integral equations. This means, in particular, that they also 
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have the same finite- volume ground state energy (as determined by the TBA approach). It is 
natural to expect that this correspondence extends to all (finite-volume) excitation energies. 

An effective way to extract information about the spectrum of the transfer matrix in lattice 
models is to use the functional relations. Using the standard fusion procedure, it is not difficult 
to show by explicit calculations^ that in the case of the dilute-A2 model considered here, the 
transfer matrix ( p.3| ) (and hence all of its eigenvalues) satisfies the relation 

T(M+7r/12)T(M-7r/12) = a{u)T [u+n / 2) + (-1)^ b{u)T{u) (4.1) 

where N is the number of sites per row (which hereafter is assumed to be even), and the scalar 
coefficients read 

a{u) = {h{u- 7r/12) h{u- 7i/Q)f , b{u) = {h{u + n/12) h{u- 7r/3)f . (4.2) 

In references Kliimper and Pearce developed a technique for transforming functional 

relations of the form ( [4.1| ) to integral equations. This is particularly useful in studying the 
eigenvalues of the transfer matrix for finite values of A^. To apply this technique one first 
requires information on the patterns of zeroes of the eigenvalues. We have used numerical 
calculations to study these patterns as follows. 

For small system sizes all solutions to the Bethe-Ansatz equations (|3.8| ) are easily found 
numerically using standard nonlinear equation solving algorithms. The eigenvalue expressions 
( p.5|) corresponding to these BA solutions, can then be compared against the results of the direct 
numerical diagonalization of the transfer matrix (|3.3|) . This enables one to numerically determine 
the patterns of zeroes of A{u) (or equivalently, patterns of zeroes of Q{u)) corresponding to each 
eigenvalue of the transfer matrix ( |3.3| ). The results obtained for the largest and a few next-to- 
largest eigenvalues of the transfer matrix T{u) are presented below. In the following, it will be 
convenient to use a new spectral variable 6, which is related to the variable u by 

e = -4i(M-7r/8-zr/2). (4.3) 

To make this change more explicit we denote the eigenvalues of the transfer matrix as T{6) = 
A(m) (assuming that 6 and u are always related by ( [4.3|) ). The eigenvalues T{6) are entire 
functions of 6 satisfying the (quasi-) periodicity relations 

T{e+Am) = T{e), T(^+4r) = exp [A^(^+4r)] T{e). (4.4) 

In fact, these eigenvalues can be written as the products 

2N 

T{e) = Ce^' l[A{i{e-e,)/A) (4.5) 

i=i 

^See also Sect. 5. 3 below. 
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Figure 1: Zeroes of the largest eigenvalue Tq{6) for = 14 and p = 0.1. 



where C and a = Yl'j=i some constants and 9j denote the zeroes of T{9). The patterns of 

these zeroes for the three largest eigenvalues (which we denote as T^^\9), T^^\9) and T^'^\9)) 
with p = 0.1, obtained numerically for = 14 sites per row are shown in Figures 1, 2 and 3a 
respectively. The zeroes 

of T(2)(^) with p = 0.3 are also shown in Figure 3b to emphasise the 
qualitative change in the position of two of the zeroes at some intermediate value of p. The total 
number of zeros in each of the figures is equal to 2N = 28. As can be seen from the figures, 
most of the zeroes accumulate near the lines 



m 9 



(4.6) 



In fact, for the largest eigenvalue T^^\9) all of the zeroes are located near these lines. The 
eigenvalue of the first excited state T^^\9) has two zeroes shifted to the positions 



9 = ±i7r + 0{p). 



(4.7) 



The eigenvalue of the second excited state T2{9) also has two shifted zeroes similar to ([4.7|), and 
in addition, has two more shifted zeroes inside the strip | '^m9\ < tt/6 which are located at 
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Figure 2: Zeroes of the 2nd largest eigenvalue Ti{6) for = 14 and p = 0.1. 

The value of 7 is controlled by the dimensionless scaling parameter 

r = 4^3 p2 N. (4.9) 

Numerical calculations show that 7 vanishes when r approaches some critical valueQ 

7 = 0, for r = re ^ 2.85 ±0.10 (4.10) 

For r > the parameter 7 is purely imaginary and approaches the value 700 = in /6 as r —>■ 00. 
For r < Tc the parameter 7 is real and increases as 0(| logr|) when r ^ 0. We found also that 
the ratio T{6)/T{6+2TTi) satisfies the condition 



T{9)/T{9+2ni) E IR, for Qm{9) = 0. 



(4.11) 



In fact, this ratio is positive for T^^\9), and for T^'^\9) when r > r^, and is neg ativefor T«(^). 
For r(2)(6') with r < r^, it is negative between the zeroes in ( [4.8| ) (z.e. for |0| < 7). Here, we 
assume that these properties of T^^\9), T^^\9) and T^'^\9) remain applicable for all arbitrarily 
large values of A^, and small p. 

^Note, that an accurate value of this critical value Vc is quite difficult to obtain due to the singular nature of 
this point which leads to severe numerical instabilities. 
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Figure 3: (a) Zeroes of the 3rd largest eigenvalue T2{6) for = 14 and p = 0.1, which corre- 
sponds to r < r^, (b) Zeroes of the 3rd largest eigenvalue T2{9) for N — 14 and p = 0.3, which 
corresponds to r > Tc. 
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For the following calculations it is convenient to pass directly to the scaling limit ( 3.3C ), 
defining the mass parameter m and the system size R as 



m 



R = N6 



(4.12) 



where S is the dimensional lattice spacing parameter. In this limit, the appropriately normalised 
eigenvalues of the transfer matrix tend to the finite values 



T{e) = lim 

N-*oo 



exp (ciA^5cosh(0/2) + C2 log cosh 0) T{9) 



(4.13) 



where the exponential multiplier in front of the limit sign is introduced for later convenience. 
The constants ci and C2 should be chosen to regularize the product over zeroes in ( [4.5|) which 
diverges in this limit. Estimating the asymptotic density of zeros 6j for large 6 from the Bethe 
Ansatz equations (|3.8|), one can show that 



1 

f 1 



Cl 



34v^ 



(4.14) 



where r = mR is the dimensionless scaling parameter previously defined in ( |4.9|) . Although it 
is also possible to calculate the value of the constant C2, this is not essential for the following 
discussion and is hence omitted here. It can be shown from the properties ( [4.4| ), that the 
functions T{9) defined in ( |4.13| ), are entire functions of the variable 9 which obey the periodicity 
relation 

T{e+Am) = T{e) (4.15) 



and have the following asymptotic behaviour for 6 ±oo 



log 
log 



T(^) 



"T Stt 



T{9)/T{e+27ri] 



~ |e±^ + 0(eT^) 



(4.16) 



where C3 is some (unknown) constant. These asymptotics are valid in the stripQ | '^m9\ < 7r/3. 
The functions T(^) are in fact the QFT counterparts of the eigenvalues of the transfer matrix 
of the lattice model. Using ( [4.1|) and (|4.13| ) it is easy to see that they satisfy the functional 
relation 



T{e+m/i)T{e-m/2,) 
Introducing the new functions 



■ cosh 6 



T{d) 



■ cosh ( 



T{d + 27rz) 



(4.17) 



Yo(^) 



cosh 6 



T(0)/T(^+27ri), Yi{e) = Yo(^+2m) 



(4.18) 



^In fact, these asymptotics are valid in a wider strip | 3m 0| < 5tt/6 — e for arbitrary small positive e. 
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one can then rewrite the functional relation ( [4.17| ) in the form 

Yo(^+m/3)Yo(e-z7r/3) 1 + Y^^O) 



(4.19) 



MO) i + Y:,\ey 

The asymptotic behaviour of the functions Yo,i(0) at ^ ±cxd is determined by ( [4. 161) to be 



log|Yo(^)| = ^e^^ + 0(e^^), log |Yi(^)| = 0(e^^), \Qme\<7i/3. (4.20) 

The derivation of the Kliimper-Pearce equations for the finite-size energies of the QFT (|2.1| ), 
is based on the following simple fact |]l| that if f{6) is a function which is regular, bounded in 
the strip | Qm{6)\ < vr/S, and satisfies the relation 

fie+tn/3) + f{e-m/3) - f{e) = g{e) (4.21) 

for some function g{9), then 

f{e)= / de' ci>{9-e')g{e') , (4.22) 



where 0(6*) is defined in ( p.9| ). 

The eqs.( ^4.21| ),( [4.22| ) combined with knowledge of the analytic properties and asymptotic 



behaviour of the functions Yo,i(^^) defined in ( [4.19|) , allows one to rewrite the functional equation 



( [4.19| ) in terms of nonlinear integral equations which generalise the TBA equations. From the 
functional equation ( [4.19|) , it is easy to see that the simplest case where Yq{6) and Yi{9) have 
no zeroes in the strip QmlO) G (—vr/S, vr/S), only requires a straightforward application of the 
lemma (see Section 4.1 below for further details). From Figures 1, 2, Sa and Sb, it can be seen 
that the only complication arising for the next few largest eigenvalues involves the function T{6) 
having a finite number of zeroes in the strip | '^'m6\ < vr/S, which are either real or occur in 
complex conjugated pairs. Suppose now, for generahty, that the function T{9) has a total of 
A + 2B zeroes in this strip where there are A real zeroes aa, a = 1, . . . ,A, and B complex 
conjugated pairs Pb ± ilb, b = 1, . . . , B. Then using ( [4. 171) and ( [4.18| ) it is easy to show that the 



eigenvalue expression ( |S.10| ) in the scaling limit ( |S.SO| ) can then be written as 

A B 

^ , ^ /vr 

lb 



E{R) = -p ^ + y e"" + 2m y cosh/??, sin - - 

ovr 2 ^ — ^ ^ — ^ V S 



a=l 6=1 

m 



2vr 



cosh^ log 1 + Yo(^)~' dO. 



(4.2S) 



In Sections 4.1-4.S below, we explicitly study the form of the resulting TBA integral equa- 
tions and energies corresponding to the ground state and first two excited state energies of the 
spectrum ( p.lO|) . 
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4.1 The TBA Equations for the Ground State 

As follows from the discussion of the properties of the eigenvalues of the transfer matrix at the 
beginning of the previous section, and as shown explicitly in Figure 1, all of the zeroes of the 
functions T{6) and T{6 + 27ri) are in this case located in the vicinity of the lines Q'm(^) = ±57r/6 
and '^m{6) = ±77r/6 respectively. In particular, they do not contain any zeroes in the strip 
I ^171(6)1 < 7r/3. Moreover, their ratio is real and positive for all real 6. Therefore the functions 

eo(^^) = logYo(e), ei(^) = logYi(^^) (4.24) 



are real analytic for | ^m{6)\ < 7r/3. Taking into account the asymptotics ( [4.20|) , it follows that 
the functions eo{6) — rcosh.6 and ei(^) are regular and bounded in the strip | '^m{6)\ < vr/S. 
Therefore one can bring [|1| the functional relation ( |4.19| ) to the form 



ej{e) = SjfimR cosh9 + ^<^j,k * l[ \e), j = 0,1 (4.25) 



fe=0 



where the kernel ^jk{0) and the function l\^\6), are respectively defined by 

= {^j,k - hk-i - kk+i) m (4.26) 

4^)(e) = log(l + e±^'=(^)), (4.27) 
the function (f){9) is defined in ( [^.9| ), and the convolution operator * is defined by 

/>oo 

f*9{e)= / f{e-&)g{e')de'. (4.28) 



Equations ( [4.25| ) are identical to the TBA integral equations ( 2.16 ). Similarly, the energy 
expression ( [4.23| ) becomes 



E,{R) = '^^^^^^-^ I cosheLl\e)de. (4.29) 

on ZTT 



which is also identical to the earlier TBA result ( p. 17 ), at least up to some appropriate bulk 



term. Consequently, the conformal properties associated with the scaling limit of these equations 
are the same as those calculated in section 2. 

4.2 The TBA Equations for the First Excited State 

As shown in Figure 2, the location of the zeroes of T{9) and T(^+27ri), is the same as above 
except that both of these functions have two extra zeroes ( [4. 81) at 6 = ±z7r. These extra zeroes 
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cancel out in ([4.18 ) resulting in the functions Yo(^^) and being real and negative for all 

real 9. Defining the functions 



eo(^) = log (-Yo(^)) , ei(^) = log (-Yi(^)) , (4.30) 

and proceeding as in the case of the ground state, one obtains from the functional relation ( [4.19|) 
and energy expression ([4.23|) , the first excited state TBA equations 

1 

ej{e) = SjflmR cosh(^) + ^ <l>j- ^ * log (l - e'^'^^^^ ) , j = 0, 1, (4.31) 

A;=0 

where the function $^^^(0) is defined in ( |4.26| ), and the corresponding first excited state energy 

rri^R log(mi?o) 



Ei{R) 



/ cosh^ log (1 - e-'°(^) ) d9. (4.32) 

2vr y_oo ^ ^ 



87r 

The leading short distance asymptotics of Ei{R) determined by these equations can be shown 



using the standard dilogarithm trick ||26| to be 

MR) ~ ^ (4.33) 

which is in agreement with the expected form Ei{R) ~ — (7r/6i?)(c — 24Ai), where Ai = is the 
second lowest conformal dimension in the operator algebra associated with the minimal CFT 
model Equations ( |4.31| ) and ( [4. 32] ) were previously conjectured in |T5[ using a different 

approach. 



4.3 The TBA Equations for the Second Excited State 

The location of the zeroes of T(^^) and T(0+27ri) in this case depends on the value of the 
parameter r. As described in the previous section (and shown in Figures 3a and 3b), most of 
the zeroes of these functions are the same as for the first excited state except that there is an 
extra pair of zeroes. The positions of these extra zeroes depends upon whether the parameter 
r < Tc or r > Tc, where the critical value Vc is specified in ( [4.10| ). For r < Tc, the extra zeroes 
of T(^) lie on the real 6 axis at the positions 6 = ±a, a > 0. For r > rc, the extra zeroes of 
T{6) lie on the imaginary 6 axis at the positions 6 = ±^7, 7 > 0. These two cases are described 
separately below. 

The case r<rc: As follows from the patterns of zeros described above, the function Yo(^) 
has two zeroes in the strip \ Q'm{6)\ < 7i/3 located on the real 6 axis at 6' = ±a, a > 0. 
Correspondingly, the function Yi(^^) has simple poles at these points. Moreover, the large 6 
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asymptotics of Yo(^) and ¥1(6*) are given by ( [4.20 ). This situation is very similar to that 
described in Reference ^ for the first excited state of the scahng Lee- Yang model. Following 
the discussion presented there, it is useful to introduce the functions 



0-0(6', a) 



tanh(| {9 -a)) 

cosh(6') — cos(?7) cosh(^) — sin(7r/ 6—1]) 
cosh(6')+cos(?7) cosh(6')+sin(7r/6— r/) 



(4.34) 



which satisfy the equations 



o-Q{9+in/3) ao{9—in/3) = 1, ai{9+in/3, rj) ai{9—m/3, rf) = ai{9, rj) 



(4.35) 



Now, defining the functions 

eo{9) = logYo(6^) - loga{9,a), e,{9) = logYi(6^) + \oga{9,a), (4.36) 

where 

a{9,a) = ao{9,a)ao{9,-a), (4.37) 

it then follows that the functions €q{9) — r cosh 6' and ei{9) are regular and bounded in the strip 
I '^m9\ < 7r/3. Applying the above lemma one then obtains from ( [4.19|) the integral equations 

eo{9) = r cosh(e) + f <p{9-9') [log (a{9, a)+e-^«(^')) - 
J —00 

- log {l+a{9, a)e-^i(^'))] d9' (4-38) 

e^{e) = rcosh^ -eo(^) 

where 0(6') is defined in ( |2.9| ). The value of a which determines the position of the zeroes of 

Yo(6^), is constrained by the functional equation (|4.17|) . Indeed, substituting 9 = a ± in/S in 
( |4.17| ) results in the condition 

Yo(a±z7r/3) = -1. (4.39) 



Using the TBA equations ( [4.38| ), and the explicit form { \i.9\ ) of the kernel 0(6'), one can rewrite 
this condition as 



'\/3sinha + 



cosh2(^-a) 
sinh2(^-a) 



log a{9, a)+e 



log l+a(9,a)e-''^^'^ 



d9' = 71 ^4/ + 2 arctan(sinh 3a)^ 



(4.40) 



where / is some integer (arising when taking the logarithm of ( [4.38| )), and where -f^^ denotes 
the principal value of the singular integral. Numerical calculations indicate that in this case 
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/ = 0. Finally, the energy expression ( [4.23 ) in the case r < becomes 
m^R log(mi?i 



E2{R) 



/oo 
cosh6' log 1 + Yo(^^)^^ 
-oo 

m'^ R \og(mRn) m f°° , ^, / n «■ \ ,^ 

/ cosh e log (a(e, a) + e"^°(^) ) dO, 

8tt 2n V / 



d9. 



Again, the leading short distance asymptotics of E2{R) can be found through an appropriate 
modification of the dilogarithm trick It is not necessary to present further details of 

this calculation here as they are well described in Appendix C of Reference . The final result 
is 



E2{R) 



97r 
lOR 



(4.42) 



which is in agreement with the expected form E2 ~ — (7r/6-R)(c— 24A2), where A2 = 1/5 is 
the third lowest conformal dimension in the operator algebra associated with the minimal CFT 
model 

The case r>rc: The consideration of this case is very similar to the case r < Tc described 
above. The function Yq{9) (and respectively Yi{9)) has two complex conjugated zeroes (poles) 
at 6 = ±i'-f, 7 > 0. Proceeding in the same manner as presented above, one obtains from ( [4.19|) 



eo(^) = r cosh(^) + logai(0,7) + 0* 
e^{e) = rcoshe-eo(^) 



Li-\9)-Ll>{9) 



(-)i 



where 



eo(^) = logYo(^), ei(^) = logYi(^). 
The value of 7 is here determined by the condition 

Yo(±^(7-7r/3)) = -1. 

The corresponding energy expression for r > is 



(4.43) 



(4.44) 



(4.45) 



E2{R) 



m?R log(mi?o 



57r 



+ 2m sin (71/3-7) / cosh ^ log ( 1 + Yo(^)"^ I dO. (4.46) 

27r V / 



5 Numerical Results and Discussion 
5.1 Numerical results 

In order to numerically determine the ground state and lowest two excited state energy levels 
Eq(R), Ei(R) and E2{R) as a function of the finite length R, it is first necessary to numerically 
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E 




r=mR 



Figure 4: The lowest three energy levels Eo{R), Ei{R) and E2{R) of the A^3,5 ± 02,i QFT 
obtained numerically from the energy expressions ( |4.29| ), ( |4.32| ), ( |4.41| ) and ( ^4.46| ), after numer- 
ically solving the respective TBA integral equations ( [4.25| ), (|4.31|) , ([4.38|) with the constraint 
i ^M ), and (lOal ) with the constraint (^M) - 



solve the respective systems of TBA integral equations for all required values of R. In the 
case of the ground state energy Eq{R), the associated TBA equations (|4.25|) are easily solved 
at particular values^] of R using a simple iterative procedure. Using these solutions it is then 
straightforward to numerically integrate the energy expression (^4.29 ) to obtain Eo{R) at these 
values of R. 

For the first excited state energy level Ei{R) defined by (|4.32| ), the associated TBA equa- 
tions ( [4.31|) can again be solved numerically at particular values of R using a simple iterative 
procedure. The only significant difference in this case, is that there are convergence problems 
which can arise due to the fact that it is now possible for the argument of the logarithmic term 
in (|4.31|) to become either very small, or even negative. Once a numerical solution of ( [4.31| ) 
is known, it is again straightforward to numerically evaluate the integral in ( [4.32|) and hence 
obtain Ei{R) for the chosen values of R. . 

For the second excited state energy level E2{R), the situation is much more complicated. 
In this case there are two distinct regions of R defined respectively by the bounds r < Vc 
^We examined a number of different R values in the range 10~^ < i? < 10. 
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First Energy Gap El -EG (TBA) 
(TCSA) 
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(TCSA) 



r=mR 



Figure 5: The two lowest energy level differences Ei{R) — Eq{R) and E2{R) — Ei{R) obtained 
numerically from the TBA data shown in Figure 6.4 and presented respectively in Table 6.1 and 
6.2, and from the TCSA method described in Section 5.4.1. 



and r > r^, where the critical value is specified by ( [4.10[ ). In fact, these two regions are 
described by quite different energy expressions and TBA equations. For r < Vc, the energy 
level E2{R) is determined by the expression (|4.41| ), which depends on the TBA equations ( [4.38|) 
and the associated constraint condition ( [4.40| ). Similarly, for r > Vc the energy level E2{R) 
is determined by the expression ( [4.46 ), which depends on the TBA equations ( [4.43 ) and the 
associated constraint condition ( [4.451) . In both of these cases, accounting for the additional 
constraint equations as well as the convergence problems (which are similar in nature to those 
described above) present some extra difficulties for the numerical solution of the two sets of TBA 
equations. However, a more significant problem in this case is that all of the aforementioned 
numerical strategies breakdown in the region around the critical value Vc of the scaling parameter 
r. This problem is avoided here through the use of an interpolating function to describe the 
energy level E2{R) in the vicinity of the crossover between these two descriptions. 

In Figure 4, the energy levels Eo{R), Ei{R) and E2{R) derived from the TBA methods 
outlined above, are plotted as functions of the scaling length r = mR. A comparison of these 
excited state TBA results can also be made with the known TCSA results [0]. In particular, it 
is most useful to compare the energy level differences Ei{R)—E(){R) and E2{R)—Eq{R) obtained 
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Table 1: A comparison of some numerical results for the lowest energy gap Ei{R)—Eq{R) 
evaluated from the TEA approach using equations ( |4.25| ), ( [4.29| ), ( [4.31| ) and ( [4.32| ), and from 



the TCSA method |JT3[ described in Section 5.4.1. 



from both of these approaches. These energy level differences are plotted in Figure 5. Some 
particular values of these energy level differences are also presented in Tables 1 and 2. It is quite 
obvious from Figure 5 that the agreement between the TEA and TCSA results is extremely 
good for the first energy level difference, but is somewhat poorer for the second energy level 
difference (especially at intermediate values of r). A total agreement with TCSA is not really 
anticipated since the dimension of the perturbation in ( |2.1| ) A = 3/4 exceeds the value A = 1/2 
above which the TCSA method requires modifications due to divergence problems. We refer the 
reader to ||3^ where these problems of the TCSA method are more thoroughly discussed. 
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Table 2: A comparison of some numerical results for the second lowest energy gap E2{R)—Eq{R) 
evaluated from the TEA approach using equations (^:25|) , (g^), ( ^381) , dOOp , (^:iTD , (|4^ ), 
( |4.45| ) and (^4.46|) , and from the TCSA method |jl3|] described in Section 5.4.1. 



5.2 Massless case 

Consider now the asymptotic expansions of the operator T(^) at large values of 6. For simplicity 
restrict ourselves to the CFT case by replacing the r cosh(6') term in ( |4.17| ), ([4.18| ) by |re^. 
The equation ( |4.17| ) can now be written as 

_r^eT{e + in/3)T{9 -in/S) 



T{9) 



(5.47) 



The ground state eigenvalue T^^\9) of T{9) does not have zeroes in the strip | Qm9\ < ii/S. 
Therefore, with an account of the asymptotics (|4.16|) , the equation (|5.4?| ) implies 



logT(°)(0) 



87r 



3'' + Ce^ + (/)*log(l + e-^°(^): 



(5.4J 



where C is an arbitrary constant and eo(^) determined by the "massless" version of (^.25| ) with 
the r cosh(6') term replaced by |re^. 
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Expanding the kernel (p.9|) in a series 



0(0) = _ 2 y gin 2TT{n+l) ^(,_2n). _ (5 49) 

TT 3 
n=l 

one gets from ( ^.48 ) 

\ogT^^\9) = -^9e' + Ce' -- V sin + Xne^^-^")' ■ (5.50) 



where 



StT TT 



Xn= I e(2"-W(l + e-^«(^))c/^ (5.51) 

oo 



The numerical values of these coefficients can be compared with the the exact results of 



where CFT's with extended W3 symmetry were studied. To facilitate this comparison let us 
introduce a new variable 

3re^ 



t' = ^ (5.52) 
where a is a normalisation constant and rewrite the first terms of the expansion ( |5.50|) as 



T(»)(t) = --|(ai)Mogi + Ct= - ^ (.i)-^/r + + . . . (5.53) 

According to [jTT| the quantities and I^^^ are the vacuum eigenvalues of the local integrals 
of motion in a particular highest weight module of the W3 algebra with the value of the central 
charge c = 6/5 (see [|n|] for further details). Their exact values [|TT| 

/f^ = --, /f = -^ (5.54) 
^ 20' ^ 10500 ^ ^ 

are in good agreement with numerical values obtained from ( p.51| ) 

/f) = -0.499999999599 10"\ 1^ = -0.115238095240 10"^ (5.55) 
Similar calculations for the first excited state of Sect. 4.2 in the "massless" limit give 

= 0.500000000002 10"\ = 0.497619047634 10^^ (5.56) 



which again in a good agreement with the corresponding exact values of |11 



^ 20' ^ 42000 ^ ^ 
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5.3 Functional Equations 

Finally note, that the functional equation ( [4.1|) can be obtained from the general set of of A^2^ 
fusion relations [32| 

Tn{u + X-'^ T„(m-A + |) = Tn-i{u) Tn+i{u) + fn+i(u-^^ T„,(m + ^) (5.58) 
where 

r_iH^O, To = (-lf /o(«-y), Ti(n)^T(n), (5.59) 

fn{u) = (p-^di{u-Vn)di{u + Vn)) , = -{{2n + l)g + 1) (5.60) 
where the theta function 9i{u) = 9i{u,p) defined in ( |A.1| ) and 



2A 

9 = 1-— (5.61) 

TT 

We would like to stress that the relation (|5.58| ) holds for arbitrary values of the parameters A 
and u in (|3.5|) , where it can be considered just as a definition the "fused" transfer matrices 
Tn{u), n > 2, in terms of Ti{u) = T{u). The eigenvalue expression ( |3.5|) and the Bethe Ansatz 
equations ( ^.8[ ) ensure that all the higher transfer matrices T„(u), n > 2, are entire functions of 
the variable u as well as T{u). For some special values of of the parameters A and u the infinite 
set of relation (|5.58| ) truncates and becomes a system of functional equations for a finite number 
of the transfer matrices T„('u). Let us a introduce a new variable 

t = p-ie2^("-3A/2) (5,62) 

and rescale the transfer matrices as 

Tn{u) = (pt)-^T„(t). (5.63) 

The relation ( p.58| ) then becomes 

T^{tq) T,(tg-i) = T„_i(t) T„+i(t) + 0„+i(t) T,(-t) (5.64) 



with 



and 



q = e'^'s (5.65) 



Ut) = [0(ztg"+5)0(-,tg— 1)]^, 0(t) = ^ (-l)V^r (5.66) 

A:=— oo 

Now it is easy to see that if the the scalar factor 0„+i(t) is omitted the relation ( [5.64| ) coincide 
with a particular case (corresponding to the A'^^ reduction) of the more general A^2^ fusion 
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relations considered in ref. |jTI[ (see eqs.(5.14), (6.47), (6.48) therein). Note, that all reasonings 
of [0 can me modified to take into account the above scalar factor 0„+i(t) in ( ^.64|) . Then 



using the determinant expressions (5.1), (5.2) of ref. (see also the discussion in Sect. 8.1 of 
that paper) one can show that if 

for some integer m then 

T„(t) = 0, Tm-k-i{t) = T,(-g'"+2t), A; = 0, 1, . . . , m - 1. (5.68) 
The equation (|4.1| ) corresponds to a particular case of this reduction 

m = 4, A = 7r/12, q = e^^'^^ (5.69) 

Another simple case is 

m = 4, A = 7r/6, g = e2^/^ (5.70) 
where ( |5.64| ) and (|5.68|) imply (assuming even) 

T(tg^) T(tg-^) = (0o(-t) + M't)) T(t) (5.71) 

The leading p — singularity of the free energy of the dilute A2 lattice model of Sect. 3 with this 
value of A = 7r/6 is log k{u) ~ logj9sin(2M), which determines to the value of the correlation 
length exponent u = 1. Therefore the scaling limit of the model corresponds to 

N ^00, p ^ 0, r = AV3pN = finite (5.72) 

where the equation ( |5.71| ) becomes 

T(^ + y) t(0-^) =2coshgcosh(0)) m (5.73) 

with 6 = logt. The massless version of this equation is related to a certain c = —2 CFT with 
extended W3 symmetry |ll| (Sect. 6. 3. 2 therein). 
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Appendix A 

The Boltzmann weights of the off-critical dilute-^M models are given by [p!8l , [27|l 

^9i(6A-m)^9i(3A+m) 



W 





a 










\ a 


a 





79i(6A)^9i(3A) 

/^(a+l) i?4(2aA-5A) S{a-l) ^i{2a\+b\)\^i{u)di{;i\-u) 
~ S{a) di{2a\+\) ^ S{a) di{2aX-\) J ?9i(6A)^9i(3A) 



W\ 



W\ 



w\ 



a ± 1 a 
a a 

a a 
a ± 1 a 

a a ± 1 
a a ± 1 



W\ 



w\ 



a a 



a a ± 1 

a a ± 1 
a a 

a ± 1 a ± 1 



u \ = p- 



^^i{'i\-u)^^{±2a\+\-u) 
^9i(3A)^94(±2aA+A) 



5(a±l)y/^ ^9i(M)794(±2aA-2A+M) 
" ' ~ ^ ' 5(a) J i9i(3A)i94(±2aA+A) 



a a 



u 



P 



?94(±2aA+3A)?94(±2aA-A)y/^ ^9i(u)^9i(3A-n 



W\ 



w\ 



w\ 



a ± 1 a 

a a =F 1 

a a =F 1 
a ± 1 a 

a a ± 1 
a ± 1 a 



?9i(2A)^9i(3A) 



^9i(2A)i9i(3A) 

5(a-l)5(a+l) i9i(u)i9i(A-u) 



52(a) 



^9i(2A)?9i(3A) 



M3X-u)M±AaX+2X+u) S(a±l) Mu)M±4aX-X+u) 

P Q /o^^ / I . ^ , o^^ + P" 



= P- 
+ P 



^9i(3A)79i(±4aA+2A) 
(3 A+m)^9i (±4aA-4A+ii) 



5(a) 'i9i(3A)'i9i(±4aA+2A) 



'(?i(3A)i?i(±4aA-4A) 

5(aTl)^9i(4A) ^94(±2aA-5A)\ i^i{u)^i{±4aX-X+u) 



S{a) 'i9i(2A) ^94(±2aA+A) J i9i(3A)^?i(±4aA-4A) 
a ^i(4aA) 



' ' ?94(2aA) 
p = /i(2A)/i(3A). 

The standard elliptic functions of variable u and nome p = e~'^ are defined as 

oo 

Mu) = Mu;p) = 2 5^(-l)V'"+'^'^^sin[(2n + 1)m] 

n=0 



^^{u) = ^i{u;p) = 1 + 2 ^(-l)V' cos(2 



nu). 



(A.l) 
(A.2) 



n=l 
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These have the useful (quasi-) periodicity properties 

i^i{u+7r;p) = -i^i{u;p), i^i{u+ir;p) = ~p~'^e~'^''''0i{u;p) (A.3) 

and 

i?4(ti+7r,p) = '(?4Kp), Mu+ir,p) = -p"^e-^^"'(?4(ti,p). (A.4) 
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